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Abstract 

The thermodynamic and kinetic properties of hydrogen adatoms on graphene are important to the mate- 
rials and devices based on hydrogenated graphene. Hydrogen dimers on graphene with coverages varying 
from 0.040 to 0. 1 1 1 ML (1.0 ML = 3.8 x 10 15 cm~ 2 ) were considered in this report. The thermodynamic and 
kinetic properties of H, D and T dimers were studied by ab initio simulations. The vibrational zero-point 
energy corrections were found to be not negligible in kinetics, varying from 0.038 (0.028, 0.017) to 0.257 
(0.187, 0.157) eV for H (D, T) dimers. The isotope effect exhibits as that the kinetic mobility of a hydro- 
gen dimer decreases with increasing the hydrogen mass. The simulated thermal desorption spectra with 
the heating rate a = 1.0 K/s were quite close to experimental measurements. The effect of the interaction 
between hydrogen dimers on their thermodynamic and kinetic properties were analyzed in detail. 

PACS numbers: 68.65.Pq, 67.63. -r, 68.43.Bc 
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I. INTRODUCTION 



The thermodynamic and kinetic properties of hydrogen (isotopes) on graphene can help us un- 
derstand the interaction between hydrogen and graphitic materials in outer space l ^ and fusion 
devices.^ Such properties also determine the realizability of the graphene(graphite)-based hy- 
drogen storage 6 and the graphene-based electronic devices.-^— And the hydrogenated graphene 
(graphite) has been one focus of scientific researches in recent years. 

Thermal desorption (TD) spectroscopy, which is also called temperature-programmed desorp- 
tion (TPD) spectroscopy, has been used in experiments to study the kinetic properties of hydrogen 
on the graphite surface. 2 ^— T Zecho et alr&H found that the saturation coverages (& sat ) of H 
and D on graphite surface are about 0.3 and 0.4 ML, respectively. The TD spectra of H and D 
(heating rate a = 1 .0 K/s) both exhibit two-peak shapes when the coverage < 0.3 x Q sat , and 
the increase in just results in some minor modification of the two-peak shapes. L. Horkekaer et 
al. found that there are mainly two kinds of hydrogen dimers responsible for the two-peak shapes 
of the TD spectra at low hydrogen coverage, which are the ortho-dimer and para-dimer. In the 
two-peak shaped TD spectra, the two desorption peaks are at 445 (490) K and 560 (580) K for H 
(D) adatoms, respectively,-!^ where the positions of the desorption peaks increase with the hy- 
drogen mass. F. Dumont et al. have simulated the TD spectra of hydrogen dimers on graphene by 
using the kinetic Monte Carlo method, 19 and reproduced the two-peak shaped TD spectra. And E. 
Gavardi et al. also have simulated the two-peak shaped TD spectra by the same method but taken 
hydrogen trimers and tetramers into consideration. 20 These experimental and theoretical results 
indicate that the hydrogen ortho-dimer and para-dimer should be still the most important dimers 
responsible for the near-two-peak shapes of the TD spectra at high coverages. This is consistent 
with that the configurations of the ortho-dimer and para-dimer always tend to be preserved in 
hydrogen trimers and tetramers. 2 ^- 

However, the kinetic properties of hydrogen dimers still need more accurate simulations, espe- 
cially the vibrations should be taken into consideration. Because the mass of hydrogen isotopes 
are so small that the vibrational zero-point energy corrections may not be negligible in kinetics. 
The thermodynamic and kinetic properties of hydrogen dimers are expected to be sensitively de- 
pendent on their microscopic structures, on which systematic and accurate investigations are still 
lacking. 

In this report, the thermodynamic and kinetic properties of H, D and T dimers on graphene 
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are simulated with a composite method consisting of density functional theory (DFT), density 
functional perturbation theory (DFPT), 23 harmonic transition state theory (hTST), 24 ~ 28 and reaction 
rate theory. 26 The effect of the interaction between hydrogen dimers on the thermodynamic and 
kinetic properties of hydrogen dimers on graphene are studied. The isotope effects in the kinetic 
properties are observed due to the inclusion of the vibrations in the simulations. The simulated 
TD spectra of H and D dimers with the heating rate a = 1.0 K/s are very close to experimental 
measurements, and the TD spectra of T (radioactive) dimers are also predicted. 



II. METHODOLOGY 



The adsorption energy of a hydrogen dimer is defined to be the energy difference before and 
after the adsorption of two hydrogen atoms on graphene, which is expressed as 



-'ads 



2E H + E GL - E GL+2H (1) 



where E H , E GL and E GL+2H are the total electronic energies of an isolated hydrogen atom, an iso- 
lated graphene layer and a graphene layer with a chemisorbed hydrogen dimer on it, respectively. 

The frequency (v) of an over-barrier jump from one local minimum state (initial state) to an- 
other local minimum state (final state) can be calculated using the quantum-mechanically modified 
harmonic transition state theory (hTST), which is expressed as 24 28 
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where v* m is the quantum-modified exponential prefactor; E ac is the activation energy; N is the 
number of atoms; co\ and a>f are the frequencies of the zth vibrational mode in the initial and 
saddle-point states in the reaction path, respectively. The numbers 3N and 3N - 1 denote that the 
initial state has 3iV real vibrational modes and the saddle-point state has 3iV - 1 real modes with 
1 imaginary mode excluded in the calculation. The activation energy here is defined to be the 
vibrational zero-point energy corrected potential barrier, which is expressed as 
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where AE P is the potential barrier in the reaction path; AF v ^(0) is the vibrational zero-point energy 
correction. When the temperature approaches to be infinite, the classical limit of the prefactor is 
expressed as 

3N 

v « = Yn^-rz (4) 

n 
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This classical-limit form is also the Vineyard's form, 25 where all of the vibrational modes are 
assumed to be completely thermo-activated. 

In this report, a mono-layer graphene is used as a structural model, which is also a safe model 
for graphite surface because the weak Van de Waals interaction between neighboring graphene lay- 
ers has negligible effect on the chemisorption properties of hydrogen . 2 - 3 ' 29 ' 30 The vacuum between 
two neighboring layers is 10 A. And 0s of 0.1 1 1, 0.063 and 0.040 ML (1.0 ML = 3.8 x 10 15 cnT 2 ) 
are computationally achieved by chemisorbing a hydrogen dimer at the center of periodic graphene 
supercells with sizes of 3 x 3 (S3), 4x4 (S4), and 5x5 (S5) times of the graphene unit cell, re- 
spectively (Fig. [T]). Besides, a para-dimer decorated 6x6 supercell (S6-P) and an ortho-dimer 
decorated 6x6 supercell (S6-0), as shown in Fig. [H are designed to allow a hydrogen dimer 
chemisorbed at the supercell center (0 = 0.056 ML). The choice of the para- and ortho-dimer 
for decoration is due to that the two dimers were found to be the most stable hydrogen dimers 
on graphene. 30 In experimental observations, hydrogen adatoms (including dimers) tend to cluster 
on graphite surface,-! 3 - making some spots covered by hydrogen while some not. This un-uniform 
distribution of hydrogen indicates that the in simulation should correspond to the effective hy- 
drogen coverage of the hydrogen-covering spots. There are 5 configurations of hydrogen dimers 
on graphene considered here, which are ortho-dimer (O), meta-dimer (M), para-dimer (P), A- 
dimer (A) and B-dimer (B) as shown in Fig. [D as well as the desorbed hydrogen molecule (H 2 ). 
Other more extended hydrogen dimers 30 are not considered. 

The structures and potential barriers are calculated using DFT and the vibrational frequencies 
using DFPT The DFT and DFPT calculations are carried out using the Quantum Espresso code 
package, 31 in which the ultrasoft 32 pseudopotentials with the BLYP 3 - 3 -^ exchange-correlation func- 
tional are used. The energy cutoffs for the wave function and charge density are 35 and 350 Ry, 
respectively. Uniform k-point grids are chosen to be 6 x 6, 5 x 5, 4 x 4 and 3 x 3 for the S3, S4, 
S5 and S6-P(0) cases, respectively. The Methfessel-Paxton smearing technique 35 with an energy 
width of 0.03 Ry is employed to speed the convergence of the numerical calculations. The elec- 



tronic density of states (DOS) for the chemisorption states in the S5 case is calculated with a 6 x 6 
k-point grid and the tetrahedron interpolation scheme. 36 The reaction paths are described by the 
minimum energy paths (MEPs) between two local minimum states, which are calculated using the 
climbing-image nudged elastic band method. 37 For the calculation of the vibrational frequencies, 
only the Y point at the Brillouin zone center is chosen. 

III. RESULTS AND DISCUSSION 
A. DFT and hTST calculations 

The calculated adsorption energies for the P-dimer (E p ds ) and O-dimer (E°. ) on the S3, S4, 
S5, S6-P and S6-0 graphene supercells are listed in Tab. HI In the S3, S4 and S5 cases, E ads 
generally increases with 0, with the only exception that the E°, in the S3 case is a little (9 meV) 
less than that in the S4 case. This indicates the attractive interaction between the neighboring 
dimers, which is consistent with the clustering of hydrogen adatoms on graphite surface observed 
in experiments M- This is because one chemisorbed dimer can cause the carbon lattice around 
wrinkled, and partly destroy the sp 2 orbital hybridization of the carbon atoms nearby, which can 
increase the affinity of those carbon atoms to bond with another hydrogen dimer. However, the E p , 
in the S6-P(0) case is larger than those both in the S4 and S5 cases, and the E° ds in the S6-P(0) 
case is larger than those in all the S3, S4 and S5 cases, although the of the S6-P(0) case is lower 
than those of the S3 and S4 cases. This is due to the effect of the periodic boundary conditions 
(PBC) we take for the DFT calculations here. In the S3, S4 and S5 cases, two interacting dimers 
are the periodic images of each other. While in the S6-P and S6-0 cases, a dimer can interact with 
the decoration dimer at the corner of the same supercell. Thus, in the S3, S4 and S5 cases, the 
PBC will partly suppress the wrinkling of the carbon lattice caused by the adsorption of a dimer, 
while in the S6-P and S6-0 cases, a dimer feels more wrinkling caused by the decoration dimer 
in the same supercell. Thus, this PBC effect on the thermodynamic stability of hydrogen dimers 
on graphene also approves that the wrinkling of the carbon lattice will result in the attractive 
interaction between dimers on graphene. 

The calculated potential barriers (AE p ) for various transitions of the hydrogen dimers on the 
S3, S4, S5, S6-P and S6-0 supercells are listed in Tab. [Ill It can be seen that AE p is dependent 
on both and the structural configuration, the same as E p , and E° ads described in the previous 
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paragraph. In all cases, the optimal reaction path for an O-dimer to be desorbed is 0-M-P-H 2 , and 
that for the P-dimer is P-H 2 , which is consistent with Hornekaer's results. 2 30 Other transitions out 
of the 0-M-P-H 2 reaction path should overcome higher potential barriers, e.g. AE p (M-B) is 0. 1 89 
eV larger than AE p (M-P) in the S6-P case (much larger in other four cases). This magnitude of 
difference will result in that the possibility of the M-P transition is about 3 orders larger than that 
of the M-B transition at 300 K (estimated by Equ. [2]). Thus, the M-B transition could be safely 
neglected when compared with the M-P transition. In addition, according to the calculations by 
Hornekaer and Sljivancanin, 23 ' 30 the evaporation of hydrogen monomers from dimers also needs 
overcome much higher potential barriers. Therefore, only the transitions in the optimal path O- 
M-P-H 2 are considered in the simulation of the TD spectra in the following. For an isolated dimer 
on graphene, there may be more than one paths that are equivalent for one kind of transition, 
e.g. there are 4 equivalent paths for the O-M transition. The number of equivalent paths for a 
transition is defined as the degeneracy of the reaction path (g pat h)- The values of the g pat hS for the 
transitions of an isolated dimer on graphene are listed in Tab. UH where the g pat hS for the adsorbing 
transitions (H 2 -M, H 2 -P and H 2 -0) are absent, because they are hardly and unnecessarily defined 
in the simulations here. However, when is high enough, the interaction between dimers will 
tend to make two different paths of the same kind transition have different AE p s, namely the 
gpath is reduced. Among the transitions in the optimal path, the P-H 2 and O-M transitions are 
the most important ones (the rate-limiting steps) for the TD spectra, which will be shown in the 
following. The A£ , p (0-M) and A£ , p (P-H 2 ) generally increase with in the S3, S4 and S5 cases, 
which indicates that it is harder for hydrogen dimers to be desorbed from graphene at higher 0s. 
This is the same as the effect of on the thermodynamic stability of hydrogen dimers on graphene, 
as discussed in the previous paragraph. However, the A£ , p (0-M) in the S6-0 case and the AE p (P- 
H 2 ) in the S6-0(P) case are larger than those in the S3, S4 and S5 cases. This is because the PBC 
effect tend to make the dimers on a S6-P(0) supercell more stable than those on the S3, S4 and S5 
supercells. Thus, it costs more energy for a dimer in the S6-P(0) case to escape from one state to 
transit to another state. 

The adsorption of a hydrogen dimer onto two C atoms in graphene will change these C atoms 
from sp 1 to sp 3 hybridized. Seen from the DOS spectra for the systems of S5+0-dimer, S5+M- 
dimer and S5+P-dimer (Fig. |2]), the hybridization between hydrogen dimers and graphene makes 
graphene be semiconducting. The band gap are 0.43, 0.36 and 0.38 eV for S5+0-dimer, S5+M- 
dimer and S5+P-dimer, respectively. However, the DOS spectra for S5+M-dimer are different 



6 



from those for S5+0-dimer and S5+P-dimer. Only in the DOS spectra of S5+M-dimer, there 
are four narrow peaks around the Fermi level, two spin-up (occupied) and two spin-down (unoc- 
cupied). And S5+M-dimer is magnetic with a spin moment of 2.0 fig, while the other two are 
nonmagnetic. The two C atoms bonded with M-dimer are equivalent in graphene, a bipartite sys- 
tem, while those C atoms bonded with 0(P)-dimer are inequivalent. The difference in magnetism 
is due to this difference in the hydrogen-bonded sites, which can be understood from the Lieb's 
theorem for the bipartite system. 38 The adsorption of a M-dimer leaves two unsaturated p z orbitals, 
which then form two quasilocal states around the M-dimer, with each state holding one electron 
and a spin moment of 1.0 fi B . This is the same as the adsorption of the hydrogen monomer on 
graphene . 29 ' 39i4 ° The contributions of the H(l s) orbitals to these quasilocal states are very small. 
However, the two unsaturated p z orbitals in S5+0-dimer or S5+P-dimer bond with each other, 
thus, there are no quasilocal states and magnetic moment in these two systems. This bonding 
between these two p z orbitals makes the O-dimer and P-dimer be more stable than the M-dimer 
(Tab. HI]). Although the transition of a hydrogen dimer changes the magnetism of the hydrogenated 
graphene, the system is dominated by nonmagnetic states. Because the M-dimer is metastable 
(seen from the AE p s in Tab. [II]), and its lifetime at 300 K is less than 1.0 s, as simulated below. 
The two hydrogen atoms of a dimer prefer to adsorb on two inequivalent C atoms in graphene. 
Likewise, from the energetic calculations for more extended dimers 30 and trimers 21 on graphene, a 
configuration is also much more stable when two nearest-neighboring hydrogen atoms are bonded 
with inequivalent C atoms in graphene. 

The calculated AF v ^(0)s for the transitions of H, D and T dimers in the optimal path 0-M-P-H 2 
are listed in Tab. [Till where the results for the S3 and S5 cases are compared. The differences 
between these two groups are very small, especially for the AF vjb (0)s of the two rate-limiting 
transitions, O-M and P-H 2 (within 6 meV). This implies that the vibrational properties of the 
hydrogen dimers on graphene are mainly determined by the localized vibrational modes, which 
is the same as those of the hydrogen monomer on graphene. 29 It also can be seen that AF vi i,(0) 
decreases with increasing the hydrogen mass. This isotope effect is also the same as that of the 
hydrogen monomer on graphene in Ref. b ! where the relationship between the isotope effects 
in phonon spectra and AF vib (0) was analyzed in detail. The AF vih (0)s for the diffusing transitions 
(O-M, M-O, M-P, P-M) are less than those for the desorbing transition (P-H 2 ), which is also the 
same as that of the hydrogen monomer on graphene. 29 The values of the AF V! j(0)s for H (D, T) 
dimer vary from 0.038 (0.028, 0.017) to 0.257 (0.187, 0.157) eV, which are large enough that can 
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not be omitted in kinetics. 

The small difference between the vibrational properties of hydrogen dimers in the S3 and S5 
cases has also been found by comparing the prefactors for the O-M and P-H 2 transitions of H, D 
and T dimers in these two cases (Fig. [3]). In the plotted temperature range, the v*. differences be- 
tween these two cases are within 20%, and the v* qm differences within 15%. Since v is just linearly 
dependent on v* m (Equ. [2]), a small difference (< 15%) in the value of v* m will not result in any 
significant deviation in the kinetic properties, e. g. the TD spectra. Consequently, the vibrational 
frequencies from the S3 case are used in the simulations of the kinetic properties of hydrogen 
dimers in all cases, which is a quite accurate treatment and brings much convenience to the nu- 
merical simulations. Seen from Fig. [31 v*. decreases with increasing the hydrogen mass, while the 
variation of v* qm with the hydrogen mass is different. The v* qm for the O-M transition is nearly in- 
variant with the hydrogen mass in the plotted temperature range, while that for the P-H 2 transition 
increases with the hydrogen mass. These isotope effects in v* d and v* qm of the hydrogen dimer are 
the same as those for the desorption and diffusion of the hydrogen monomer on graphene, 29 which 
has been understood from the isotope effect of the spectra of the localized vibrational modes of 
hydrogen and analyzed in detail in Ref. The hydrogen dimers are kinetically active in the 
temperature range from 350 to 650 K, above which the dimers have been completely desorbed 
(shown below). In this temperature range, the v* m s for the O-M transition vary from 8 to 15 x 10 12 
s" 1 , and those for the P-H 2 transition vary from 5 to 9 x 10 12 s" 1 . 

The calculated jump frequencies for the transitions of H, D and T dimers on the S6-P supercell 
are shown in Fig. |4](a - e). The jump frequency of a transition in other cases (v x , X = S3, S4, S5, 
S6-0) can be deduced from the corresponding v in the S6-P case by the equation 

v* = v^exp(-/l) (5) 
k B T 

where 8e is the difference of the AE P for a transition in the X case and the S6-P case. The values 
of the AE p s for all the transitions in all the cases are listed in Tab. [Ill The isotope effect in v is 
that it decreases with increasing the hydrogen mass, which is due to the isotope effect in AF V!fo (0) 
(Tab. InTl) and is the same as that of the hydrogen monomer on graphene 29 . In the temperature 
range from 350 to 650 K, the jump frequencies for the M-0 and M-P transitions are at least 
4 orders larger than those for the O-M, P-M and P-H 2 transitions. Because v is exponentially 
dependent on E ac or AE P (Equ. |2] and H]) and the AE p s for the M-0 and M-P transitions are 
much less than those for the O-M, P-M and P-H 2 transitions. The life time of the M-dimer (t m ~ 
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vm p+v m ^ e case is less than 1 .0 s at 300 K, and decreases exponentially with temperature. 
And tm is much less in other cases due to their lower AE p s for the M-P and M-0 transitions 
(Tab. HU). Thus, although the M-dimer is an important configuration in the kinetic transitions 
of a hydrogen dimer on graphene, it has not been observed in experiments^ and theoretical 
simulations.— Furthermore, Vp_# 2 transition is at lest 2 orders larger than vo-m, which indicates 
that the O-dimer is kinetically more stable than the P-dimer on graphene. 



B. The simulation of the TD spectra 

In the optimal path of 0-M-P-H 2 , the jump frequencies for the M-0 (y M -o) an d M-P (y M -p) 
transitions are at least 4 orders of magnitude larger than those for other transitions in the temper- 
ature range from 350 to 650 K. Thus, it is reasonable to assume that when a dimer succeeds in 
jumping from O-dimer state to M-dimer state, it will immediately turn to be O-dimer or P-dimer 
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for the diffusion of lithium in the interca- 



state. Following the treatment by Toyoura in Ref. 
lated graphite, the mean time of the total jump is the sum of those of the individual jump steps, 
r = Ti + r 2 , where t x = (g -MV -M)~ l and r 2 = (g M -pV M -p + gM-d»M-o)~ x , with the g path of 
each transition path being accounted. An O-P transition is accomplished if the M-dimer finally 
succeeds in jumping to be a P-dimer. Then the jump frequency for this O-P transition is given by 

gM-pVM-P 1 

V -P = X (6) 

gM-OVM-O + gM-pVM-P T\ + T 2 

VM-P ,„s 

= gO-P VO-M (V) 

VM-0 + Vm-P 



where Equ. [7]is exact only for the isolated dimers on graphene and g _ P (= go-MgM-p/2) equals 
4 there; sm-pvm-p Qr — vm__p — - ls ^ e success probability of the M-P jump from a M-dimer. 

gM-O v M-0+gM-pVM-P Vm-O+VM-P r J 

g -p is regarded as the path degeneracy of the O-P transition, and Equ. [7]is used for non-isolated 
dimers on graphene where go-p is treated as a variable. Likewise, the jump frequency for the P-0 
transition is given by 

V M-0 fQ , 

Vp-0 — gP-0 Vp-M \o) 

Vm-0 + Vm-P 

The interaction between neighboring hydrogen dimers tends to lower the values of g ^ P and g P _ 
from 4. Both of g -p and g P _ are taken to be equal in the simulation of the TD spectra for 
convenience. In Fig. |4] (f) and (g), the calculated v -ps and v P - s for H, D and T dimers are 
shown, with g _ P and g P - both being taken an intermediate value of 2. v -p and vp_ are mainly 



determined by v -m and v P - M , respectively, with the multiplying factors in Equ. [7J and [8] just 
varying from 0.5 to 1.0 in the temperature range from 350 to 650 K. vo-p and vp_o decrease with 
increasing the hydrogen mass, the same isotope effect as other transition frequencies in the optimal 
path of 0-M-P-H 2 (Fig. ©. 

The ratio of O-dimer and P-dimer (n : n P ) at the starting moment (t=0) was reported as 
0.152, 0.222 and 0.47 19 , respectively. In this report, we take this ratio as 0.2 (n = 1.0 and n P = 
5.0). If the interaction between hydrogen dimers approximately acts as an average background 
effect, the reaction of the hydrogen dimers on graphene is first-order, which is consistent with 
experimental observations.-^*^ The reaction-rate equations for the TD spectroscopy of hydrogen 
dimers on graphene can be expressed as 

dn (t) 



dt 
dn P (t) 

dt 



= -v -p(T)n (t) + Vp-o(T)n P (t) (9) 
= -[v P -o(T) + vp„ H2 (T)]np(t) + v ~p(T)n (t) (10) 



n (t + A?) = n (t) + ^^-At (11) 

dt 

n P (t + At) = n P (t) + ^£^At (12) 

dt 

where At is a time interval. The desorption rate of hydrogen dimers, or the hydrogen-gas yield per 
unit time, is defined as 

R d es(T,t) = vp- H2 (T)n P (t) (13) 

where T = at and the heating rate (a) is set as 1.0 K/s in accordance with experiments . 2 i 16 i 1 7 

The simulated TD spectra for H, D and T dimers in the S6-P case are shown in Fig. [5] (a). 
The two-peak shaped TD spectra have been reproduced, with the desorption peaks being at 436 
(460, 469) and 562 (574, 581) K for H (D, T) dimers, respectively. These simulation results are 
very close to experimental meas urements h 1 61 1 7 The kinetic processes behind the TD spectra are 
described in the following. As the temperature increases (T= 1.00, the P-dimers are first desorbed 
through the P-H 2 transition above 350 K. At this time, there are negligible amount of P-dimers 
transiting to be O-dimers, because v P ^ M (vp-o) is a t l est 2 orders less than v P ^ Hl . Afterwards, 
the O-dimers are reduced at higher temperatures (> 475 K) through the 0-M-P-H 2 transition. 
The variations of n P (T/a) and n (T/a) are shown in Fig. [5] (b). When a dimer succeeds in 
transiting from an O-dimer through a M-dimer to a P-dimer, it stays as a P-dimer for a short 
time of about l/vp_ ff2 seconds before being desorbed, which is reflected by the small occupation 
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number of the P-dimers (n P (T/a)) at the temperatures around the 2nd desorption peak (the inset 
in Fig. [5] (b)). The mobilities of the hydrogen dimer in diffusion and desorption decrease with 
increasing the hydrogen mass, which is due to the isotope effect in v or AF vib (0). In experimental 
measurements , 2 ' 16117 the positions of the two desorption peaks in the TD spectra for H (D) dimers 
are at 445 (490) K and 560 (580) K (the average < 0.12 ML), while those in the simulated 
spectra are 436 (460) K and 562 (574) K (the effective = 0.056 ML). The simulated isotopic 
differences of the desorption peaks of H and D dimers are 24 and 12 K for the 1st and 2nd peaks, 
respectively, in comparison with the experimental measurements of 45 and 20 K, respectively. 
These little deviation between our theoretical simulations and experimental results maybe come 
from that a small amount of uncertainty has been introduced in experiments when estimating the 
desorption-peak positions from the somewhat chaotic TD spectra of H, 17 or from the isotope effect 
in the effective of the hydrogen-covered spots on graphene in experiments. The latter possibility 
is understandable since that the thermodynamic stability of the hydrogen dimer on graphene should 
be hydrogen-mass dependent, thus, the distribution of hydrogen dimers on graphene should be 
hydrogen-mass dependent. This is like the isotope effect of & sat , which are about 0.3 and 0.4 ML 
for H and D on graphene, respectively. 17 The difference in the effective for H and D dimers 
results in the difference in the interaction between dimers and then the difference in AE p . Based 
on this point, the isotope effect in the experimental TD spectra of hydrogen dimers consists of the 
contributions from the isotope effects in v and 0. The simulations here, however, just consider the 
isotope effect in v. 

In order to observe the effects of and the structural configuration on the TD spectra, the AE p s 
for the transitions of O-M, M-O, M-P, P-M and P-H 2 in all the cases considered here (Tab. HU) are 
used to simulate the TD spectra. The positions of the two desorption peaks in each case are listed 
in Tab. [IV] The 1st desorption peak shifts nearly linearly with the AE P of the P-H 2 transition by 
34 K/0.1 eV, while the 2nd peak does not have a simple monotonic variation due to the complex 
transitions of hydrogen dimer behind but still generally increases with AE P . Additionally, the 
interaction between dimers also tends to lower the degeneracy of the diffusing reaction paths, g -p 
and gp-o, which are taken equally here and together denoted as g pat h below. The position of the 1st 
desorption peak is obviously not dependent on the diffusing g pat i,. Seen from Fig. |5](c) whereas, the 
positions of the 2nd desorption peaks for H, D and T dimers are nearly logarithmically dependent 
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on gp at h, which is approximately expressed as 

T(g path ) = T(l)-36xlg(g path ) (14) 

where T(l) is the peak position when g path = 1. And T{g path ) increases by about 1 1 K with g path 
decreasing by a half. 

From the analysis above, the interaction between hydrogen dimers on graphene is coordinated 
with the effective hydrogen coverage and the structural configuration. Such interaction affects 
the potential barriers and the reaction-path degeneracies of the hydrogen-dimer transitions on 
graphene. As a result, the jump frequencies and TD spectra vary with the hydrogen coverage 
and the structural configuration. However, the localized vibrational properties of hydrogen dimers 
are not significantly influenced by the hydrogen coverage. 

IV. CONCLUSIONS 

The thermodynamic and kinetic properties of H, D and T dimers on graphene have been simu- 
lated with a composite ab initio method consisting of density functional theory, density functional 
perturbation theory, harmonic transition state theory, and reaction rate theory. The effects of the 
hydrogen coverage and the structural configuration on the thermodynamic and kinetic properties 
have been revealed by varying the coverage from 0.040 to 0.111 ML (1.0 ML = 3.8 x 10 15 cm -2 ). 
It has been found that the interaction between hydrogen dimers influences the effective hydrogen 
coverage and the configuration of hydrogen dimers deposited on graphene and affect the thermo- 
dynamic and kinetic properties (E ads , AE p , v(T) and Rd es (T, t)) of those dimers. The vibrational 
zero-point energy correction (AF vih (0)) and the jump frequency (v) both decrease with increasing 
the hydrogen mass. However, the isotope effect in the positions of the desorption peaks in the TD 
spectra is inverted. In a word, the mobility of the hydrogen dimer in desorption and diffusion de- 
creases with increasing the hydrogen mass. The simulated TD spectra (heating rate a = 1.0 K/s) 
are quite close to experimental measurements, and the effect of the dimer interaction has been 
clarified in detail via the potential barrier AE p in various cases and the reaction-path degeneracy 

gpath- 
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TABLE I: The hydrogen coverage (0) and the adsorption energy of the para-dimer (E^ ) and ortho-dimer 
(E°, ) in various cases (1.0 ML= 3.8 x 10 15 cm" 2 ). 



Supercell (ML) E^ s (eV) E° ads (c\ 

53 0.111 2.714 2.555 

54 0.063 2.562 2.564 

55 0.040 2.538 2.518 
S6-P 0.056 2.604 2.758 
S6-0 0.056 2.578 2.652 



TABLE II: The reaction-path degeneracy {g pa th) for an isolated hydrogen dimer on graphene, the potential 
barrier (AE P ) for the transitions of a hydrogen dimer on various graphene supercells. 

AE p (eV) 



Path 


Spath 


S3 


S4 


S5 


S6-P 


S6-0 


O-M 


4 


1.679 


1.699 


1.664 


1.690 


1.713 


M-0 


2 


0.405 


0.510 


0.596 


0.689 


0.609 


MP 


2 


0.351 


0.538 


0.607 


0.773 


0.660 


P-M 


4 


1.792 


1.711 


1.674 


1.620 


1.690 


P-A 


2 


1.821 


1.637 


1.602 


1.605 


1.636 


A-P 


1 


0.223 


0.449 


0.525 


0.558 


0.587 


M-B 


2 


0.605 


0.806 


0.836 


0.962 


0.912 


B-M 


2 


1.106 


1.139 


1.110 


1.121 


1.113 


0-H 2 


1 


2.280 


2.228 


2.360 


2.214 


2.262 


H 2 -0 




4.414 


4.298 


4.470 


4.191 


4.345 


M-H 2 


1 


0.801 


0.934 


0.994 


1.149 


1.039 


H 2 -M 




4.191 


4.187 


4.143 


4.126 


4.227 


P-H 2 


1 


1.456 


1.358 


1.336 


1.461 


1.414 


H 2 -P 




3.398 


3.448 


3.461 


3.593 


3.573 
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TABLE III: The vibrational zero-point energy correction (AF v $(0), in eV) for the transitions of H, D and T 
dimers in the optimal path (O-M-P-H2) in the cases of S3 and S5. 

S3 S5 



Path H D T H D T 

O-M 0.141 0.107 0.093 0.135 0.102 0.088 

M-0 0.038 0.028 0.017 0.051 0.034 0.027 

M-P 0.053 0.032 0.023 0.079 0.053 0.042 

P-M 0.138 0.106 0.093 0.144 0.111 0.075 

P-H 2 0.257 0.187 0.157 0.251 0.184 0.156 



TABLE IV: The positions (in K) of the two desorption peaks in the TD spectra of H, D and T dimers in 
various cases (a - 1.0 K/s). 

1st desorption peak 2nd desorption peak 



Supercell H D T H D T 

53 434 458 468 534 536 551 

54 401 424 434 551 563 569 

55 393 417 427 537 548 554 
S6-P 436 460 469 562 574 581 
S6-0 420 443 453 562 573 579 
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FIG. 1: (Color online) The structures of the periodic graphene supercells with the sizes of 3 x 3 (S3, 18 C 
atoms), 4x4 (S4, 32 C atoms), 5x5 (S5, 50 C atoms) and 6x6 (S6, 72 C atoms) times of the graphene unit 
cell. A hydrogen dimer is allowed to be chemisorbed at the center of each supercell (marked by a purple 
circle). There are 5 configurations of hydrogen dimers on graphene considered here, which are ortho-dimer 
(O), meta-dimer (M), para-dimer (P), A-dimer (A) and B-dimer (B), as well as the desorbed hydrogen 
molecule (H2). And the S6 graphene supercell is decorated with a para-dimer (S6-P) or an ortho-dimer 
(S6-0) at the corner (shaded in gray) of the supercell. The gray spheres are carbon atoms, and the blue 
smaller spheres are hydrogen atoms. 
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FIG. 2: (Color online) The DOS spectra of the S5 graphene supercell with an O-dimer, M-dimer and P- 
dimer. The vertical line at the Fermi level guides the eyes. 
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FIG. 3: (Color online) The variation of the v* qm and v* cI with respect to temperature in the S3 and S5 cases 
for (a) O-M and (b) P-H2 transitions of H, D and T dimers. The prefactors for the S3 case are plotted with 
thicker lines, and those for the S5 case with thinner lines. 
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FIG. 4: (Color online) The variation of v with respect to the inverse of temperature for the (a) O-M, (b) 
M-O, (c) M-P, (d) P-M, (e) P-H 2 , (f) O-P and (g) P-0 transitions of H, D and T dimers. The corresponding 
v* m s are used as the prefactors. 
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FIG. 5: (Color online) (a) The TD spectra of H, D and T dimers {a - 1.0 K/s) with the starting dimer 
ratio «o(0) : «p(0) = 1:5. (b) The variations of np(T/a) and no(T/a). The inset shows the details of the 
variation of wp(T/or) (in the logarithmic scale) from 450 K to 650 K. (c) The variations of the positions of 
the 2nd desorption peak in the TD spectra for H, D and T dimers with respect to the degeneracy (g pa th) of 
the diffusing reaction path . g pa th is plotted in the logarithmic scale. 



20 



